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Abstract 

The widths for pionic decays of lowest two excited doublets (0 + , 1 + ) and 
(l + ,2 + ) of heavy mesons are studied with QCD sum rules in the leading 
order of heavy quark effective theory. The ambiguity due to presence of two 
distinct 1 + states are solved. The soft pion approximation is not used. 
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I. INTRODUCTION 



The properties of excited heavy mesons have attracted attention in recent years. The 
mass and decay widths have been studied with potential model combining with data for K 
mesons in with relativistic Bethe-Salpeter equation in |2||| and with QCD sum rules 
in [|J. In work || the masses of low lying excited heavy mesons (2 + , 1 + ) and (1 + , + ) were 
studied with QCD sum rules in the heavy quark effective theory (HQET) || up to the 
order of 0(l/m,Q), where rriQ is the heavy quark mass. The present work is devoted to the 
decay widths of these mesons in the same formalism. 

One problem encountered in studying the decay widths of excited heavy mesons with 
QCD sum rules is the following. Except for the lowest states 0~, 1~, the spectra contains 
a pair of states for any spin-parity j p with close values in their masses but quite different 
in magnitudes of their decay widths. In our case, one of the two 1 + states is a narrow 
resonance decaying mainly by emitting a D wave pion, while the other one is a very wide 
resonance decaying by emitting a S wave pion. An interpolating current used for the 
narrow 1 + state with a small coupling to the other 1 + state may cause sizable error in the 
result of calculation. It is only in the iuq — > oo limit, there is a conserved quantum number 
ji, the angular momentum of the light component, which can be used to differentiate the 
two states. Therefore, HQET has important advantage for this purpose. 

The proper interpolating current J^pj" 3 for the states with the quantum number j, P, 
je in HQET was given in ||. They were proved to satisfy the following conditions 

(<> ./;:Kr'(0). /'./"../;) = f P ^p P >5 kj <f^ , m 

i (0|T (^r^)^p"5'( )) 1°) = S jf 6 Pp/ 6 Jlfe (-iySgr Pl ■ ■ -g?* 

x J dt6(x-vt)IL P j t (x) (2) 

in the vhq — > oo limit, where t] ai '" aj is the polarization tensor for the spin j state, v 
is the velocity of the heavy quark, = g af3 — v a v 13 is the transverse metric tensor, S 
denotes symmetrizing the indices and subtracting the trace terms separately in the sets 
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(aii ■ ■ ■ oij) and {(3\ - ■ ■ j3j), fp t j e and Hpj e are a constant and a function of x respectively 
which depend only on P and ji. Because of equations (HD and (|]), the sum rule in HQET 
for decay widths derived from a correlator containing such currents receive no contribution 
from the unwanted states with the same spin-parity as the states under consideration in 
the m,Q — > oo. Starting from the calculations in the leading order, the decay amplitudes 
for finite tuq can be calculated unambiguously order by order in the I/tuq expansion in 
HQET. 

II. SUM RULES FOR DECAY AMPLITUDES 

In the present work we shall confine ourselves to the lowest lying excited states in the 
leading order of 1/rriQ expansion. Denote the doublet (l + ,2 + ) with ji = 3/2 by (Bi, B 2 ) 
and the doublet (0 + , 1 + ) with ji = 1/2 by (B' ,B[). From covariance and conservation of 
the angular momentum of the light component in the rriQ — > oo limit, the amplitudes for 
the decay of B%, B\ to the ground states B, B* in the doublet (0~, 1~) with ji = 1/2 have 
the following forms 

m(b 1 - zrvr) = / e>(« - IgrqMBi, B*) , (3) 

M(B* 2 Bit) = I V , v q?q»g(B* 2 , B) , (4) 
M(B* 2 B*ir) = I ie aPap e* a v^q^g(B* 2 , B*) , (5) 

where r]^, rj^ and e M are polarization tensors for states 2 + , 1 + and 1~ respectively. q t ^ = 
% — v ' Q v fj.- I — V%, 1 for charged and neutral pion respectively. It can be shown M 
by combining heavy quark symmetry and chiral symmetry that the constants in (^|)-(^|) 
satisfy 

g(B* 2 , B) = g(B*, B*) = j| g(B 1 , B*) = Q . (6) 

For deriving the sum rules for the constants in the decay amplitudes @ and fl5|) 
we consider the correlators 
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J d 4 x e-*»<7r(g)|r (^(0)^(0:)) |0> 

= J Jd 4 xd 4 y e-^^^^irrjf^j^i^jJCO)^^*)^ |0) 

= (tfrf-^^GftB.Cw.w'), (7) 

/ d 4 x e-*-(7r(g)|T (v.^^O*)) |0> = (tfgf - t^\ 2 ) G b *b(^') , (8) 

| rf 4 x e- tt -(7r( ff )|T (^(0)^(2;)) |0) = % - (e^v p q ta qf + a <-> (3) G B * B *(cu,u;') , (9) 

where k = k' + q, oj = 2v ■ k, u' = 2v ■ k' and q 2 = 0. PCAC is used in the first equality of 
(0). The forms of the right hand side of (|7|), (§) and @ are determined by that a, (3 and 
7 are transverse indices, x — y = v t on the heavy quark propagator and that the pion is in 
the D wave state due to conservation of angular momentum of the light component. The 
interpolationg currents are given in || as 



J 



fa 

1-lS 



2 



/3 ~ 
4 



(10) 



< a + T =\\\ h v { ~y- + 1?V? - \g?^f t )q, (11) 



j t,\ = v \ ' j d,-4 = v ^ ^ ' ( i2 ) 

where fo„ is the heavy quark field in HQET and 7^ = 7^ — t^. 

Let us first consider the function GbxB* (w, oj') in (0). As a function of two variables, it 

has the following pole terms from double dispersion relation 

f_if + zg(BxB*) c c ' 

_ ■2' / +'2 yV _ 1 >_ l_ t I ° (1$) 

(2A_ ) i -cu / )(2A +j | -w) 2A_i-o/ 2A +) |-w' v ; 

where fpj e are constants defined in (P, Ap Jf = mp^ — tjiq. As explained in Section |, 
only one state with j p = 1 + contributes to ( jl3|) as the result of equation ([[]). This would 
not be true if the last term in ( |T0D is absent. 

Usually, the soft-pion approximation u = u' is used in QCD sum rules for decay 
amplitudes. This is not a good approximation since the difference u — u' = 2v ■ q is 
about 1 GeV in the real decay process as can be seen from the observed spectra of D 
system and the calculation results of ||. We instead put 



w-u/ = 2(A +j 3 - A_ i) = 2A . (14) 

This is satisfied by the real decay process in the leading order of 1/itlq. We then use 
dispersion relation of one variable to Gb 1 b*{ uj , the pole terms of which are 

(2A_ j i -uo'f (2A_ ) i -u') ' 1 J 

f 1 + i 

The free heavy quark propagator in HQET is J dr S(x — vr) — - — . The interaction of 
the heavy quark with the gluon field A in the leading order of 1/ttlq expansion is ghv ■ Ah. 
We shall use the Fock-Schwinger gauge x^A^i^x) = throughout this paper. In this gauge 
the full propagator of the heavy quark (0\T(h(0)h(x))\0) in the leading order of 1/mq is 
identical to the free one, because v = x/r on the heavy quark line. Therfore, the correlator 
(0) has the form 

dr e^ T tr{r(7r(g)|AgWg(0)|0)} = (fcV - 9t aP ) I dr e^ T f(vq, r), (16) 



where T is a vertex factor containing 7 matrices. Since in our approach v ■ q is fixed to 
be equal to A by ( |T4] ) and the invaraint function G(BiB*) is considered as a function of 
a single variable u (or u'), the integral in r is dominated by small r, that is, the small 
distance x = v ■ r for u in the deep euclidian region. Therefore for deriving QCD sum rules 
we can calculate the correlator (|7|) by the operator product expansion. This is an expansion 
over the small quantity r, whch after the integration becomes an assymptotic expansion 
in powers of 1/u. By dimensional reason it is a power series in v • q/u and Aqqd/^ 
which is meaningful for sufficiently large value of to. After Borel transform the series is 
convergent when v ■ q/T is sufficiently small, where T is the Borel transform parameter. 
We emphasize that this is so because v ■ q is fixed in our approach. The perturbation terms 
do not contribute since they do not violate the conservation equation d^j^ = 0. For the 
higher terms we expand the quark field q(x) in the matrix element (n(q)\T(q(x)q(0))\0) . 
The matrix elements appearing in the expansion have the general form (0), fl38p and 



in Appendix. The terms in (p7| )-(p9D which contribute to the QCD sum rules in the soft 



pion approximation were used in Q. Here we keep all terms in these formulas in the 
calculations in order to avoid the soft pion approximation. Substituting these equations in 
the correlator (|7]) and making Borel transform to supress the contribution of higher states 
we find the sum rule in the leading order of HQET. 



g{B X B*) -2A i/T , -2A i/T 

7f e -a +(c + c)e ^ + 



\/|b + |(c 2 -rf 2 -e 2 A)]e 2 ^. 



(17) 



In the above equation we have neglected all radiative corrections and operators with di- 
mension > 6. In the derivation of above equation we need to rotate the integration path 
clock- wise to the imaginary axis of r and redefine the integral variable r' — ir. The result 
is then obtained by using the Borel transform formula B^ J dTe lTUJ r n = (— i) n+l (l/T) n in 
the negative region of uj. 



The dots in the left hand side of (|17|) represents terms from the contributions from the 

radial excited states of the same j p and je and the continuum. In order to estimate these 

terms we make the duality assumption that above u = uo c the dispersion integral of the 

hadron side is equal to that obtained from the QCD calculation. The spectra function 

Pqqd(uj) of the latter is obtained by performing the Borel transform Bf^ to the right 

hand side of (|T7j). The result for Pqqj}{uj) contains two terms proportional to 5(uj — 2 A) 

r<x> pQQj^{uj')duj' 

and S'(cu — 2 A) respectively. This implies that / — vanishes for u c > 2 A ~ 1 

Jui c lo — to 

GeV. Therefore we shall neglect the contributions from the radial excited states and the 
continuum to the sum rule (|T7j). 

We can multiply ( |lTD by exp(2A_ i/T) and differentiate in 1/T in order to eliminate c 
and c'. By using (^) we obtain the following sum rule for the constant Q 

Gf-if + ,l =2{A + sg 2 + (c 2 -d 2 -e 2 A)(l + ^)}e . (18) 

The sum rules for the correlators (|8|) and ([]) can be obtained by similar calculations. 
Equations similar to (|T^) are obtained and we have checked that the relations in ( ^ ) are 
satisfied. 



The decay amplitudes of the members of the doublet (J, +,1/2), B' Q and B[, have the 
following form 

M(B' — > Btt) — I Q' , (19) 

M(B[ -> B*tt) = I e* ■ j]Q' (20) 

in the mq — > oo limit. The sum rule in HQET for the coupling constant can be derived 
similarly. From ||, the interpolating currents for B' Q , B[ can be chosen as 

J oV,| - 71 M > = ^ ^ 9 • (21) 

As pointed out in ||, for any state (j,P,jg) besides the interpolating current with the 
lowest dimension there is another one which has an additional factor —ipt inserting before 
q. For the doublet (J, +, 1/2), they are 

J' ] + ,i = KHWtQ , J'l+,i = ^ h vl 5 ^{-i)p t q . (22) 

As discussed in ||, there is some motivation for using the latter currents for the corre- 
sponding constant f' + 1 , 2 survives in the non-relativistic limit of the Bethe-Salpeter wave 
function, while the constant f+,1/2 for the current (|2"ID vanishes in this limit, since B' and 
B[ are P wave states in the non-relativistic quark model. 

After Borel transformation and eliminating the single pole terms by differentiation, the 
sum rule derived by using the interpolating current fl2"2"|) is 



/_ i/i i0' = 2[A+ 1 (-# 2 A 2 + 3h 2 A + 3hi) - (-3&i + 3 Cl A + 3diA 

+ 4c 2 A 2 + 2d 2 A 2 -e 2 A 3 )(l + ^±^)]e , (23) 



where A = A+i/ 2 — A_ 5l / 2 . If the interpolating current ( f2~T| ) is used instead of (p2[), the 
sum rule turns out to be 



f M = I (-A + , 1Q - § ) - ^ + ^£ ) , (24) 



where 



a = 2{h + / 2 A) , p = 4{ 9l A + g 2 A 2 + h) , 

7 = 4(ai + a 2 A + 2ciA + ei A 2 + e 2 A 3 ) (25) 

and A = A+,i/ 2 - A-,i/2- 

III. DETERMINATION OF THE PARAMETERS 

In order to obtain the decay constants Q and Q' from (0), (|23D and (|4]) we need to use 
the mass parameters A's and the coupling constants /'s of the corresponding interpolating 
currents as input. A ]/2 and f-,1/2 can be obtained from the results in || as A 1/2 = 0.5 
GeV and /-,i/2 — 0.24 GeV 3 / 2 at the order a s = 0. Notice that the coupling constant /_,i/2 
defined in the present work is a factor l/y/2 smaller than that defined in A +j 3/ 2 , h+,1/2 
and A' + are given in 0. f+,3/2, f+,1/2 and /+ !/ 2 can be determined from the formulas 
(34), (27) and (28) of reference f| derived from sum rules for two point correlators. The 
results are 

A +j3/2 = 0.95 GeV / +)3/2 = 0.19 ± 0.03 GeV 5/2 , 

A' +jl/2 = 0.9 GeV f' +1/2 = 0.28 ± 0.03 GeV 5/2 , 

A +jl/2 = 1.05 GeV / +) i /2 = 0.32 ± 0.03 GeV 3/2 . (26) 

(|40|)-(^5|) still contains several parameters, wig is known to be 0.8 GeV 2 . m\ was 
calculated in || with QCD sum rules with the result m 2 = 0.2 GeV 2 . We have obtained 
QCD sum rules for e\ and e 2 by considering the two point functions 

i J e^c/ 4 x(0|T{g(x)z 75 yg(x)g(0)y (n • D) 2 * 75 g(0)|0> = (n ■ q) 2 A(q 2 )5 ab , (27) 

1 J e ill x d A x(0\T{q(x)ii l5 ^q(x)q(0)j(n ■ Df^ 5 q(0)\0) = (n • q) A B{q 2 )5 ab , (28) 

where n 2 = 0. The resulting QCD sum rules are 

X ( ™ >0 = i^/o e 1 ds+ 2^ {a ° G ) ° + ir™^- (29) 



20 / m 2 A \ (qq)l , . 

In obtaining (j30"l) we have eliminated the a\ pole term by differentiation with respect to 
the Borel parameter 1/M 2 . vha — 1.3 GeV is the mass of a\. In the numerical calculations 
for ei and e 2 we use the following standard values of the vacuum condensates: 



, = -(0.24 GeVy , 
(a s GG} = 0.038 GeV 4 (31) 

as well as a s = 0.4. From ( p9|) and (|30|) we obtain the following numerical values for e\ 
and e2, 

ei = -0.014 ±0.002 GeV 2 , 

e 2 = 0.015 ± 0.002 GeV . (32) 



IV. NUMERICAL RESULTS AND DISCUSSION 



We now turn to the numerical evaluation of the sum rules for Q and Q 1 . By using (|26| 



and ([H]) together with the values of parameters determined above we obtain Q and Q' as 
functions of T from the sum rules flT8|) , ( ^3|) and (p4|). The results are plotted as curves in 
fig. 1. The lower limit of T is determined by the requirement that in the sum rules before 
differentiating with respect to 1/T (such as (|17D), the terms of higher order in 1/T in the 
operator expansion are reasonably smaller than the leading term, say < 1/3 of the latter. 
This leads to T > 1.2 GeV and T > 3 GeV for the sum rule (p~8| ) and ( |23"D respectively. 

Usually the value of T is constrained from the upper side by the requirement that the 
contribution from higher states is sufficiently suppressed. The continuum model discussed 
below equation (|T7]) allows T to be any large value. However, in order that the results 
depend less on the continuum model it is good if T is not much larger than the 2(A^ 1 - ) — A), 
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where mg + is the mass of the radial excited state. Quark model calculations [T(J 
indicate that A*^ — A is roughly 0.7 GeV. Therefore, this requirement can be satisfied for 
the sum rule for Q. We find within the stability window subject to the constrain disussed 
above 

Q = 5.7 ±0.6 ±0.7 GeV" 2 , 1.3 GeV < T < 1.5 GeV. (33) 

From the view point stated above, the sum rule (^) for Q' is not very satisfactory since 
the lower limit of T is too high. If we somewhat arbitrarily take the width of the window 



to be 0.3 GeV, we find from (23) 



Q' = 1.0 ±0.1 ±0.5, 3 GeV<T<3.3 GeV. (34) 

In (^) and (|34l) the first error refers to the variation with T in this region, while the second 
error takes into account the uncertainty in the values of the coupling constant /'s. The 
errors quoted above do not include those from higher order terms in OPE and the I/tuq 
expansion and the values of parameters characterizing the condensations. 

However,the (|33|) and (|34"D have another source of uncertainty. These sum rules depend 



on the parameter e\ which is obtained by the sum rule (|29|). It is well known || that the 
sum rules derived from the correlator of two pseudo-scalar currents give too small value for 
the physical quantities. In the case of two pseudo-scalar current without derivative, the 
phenomenogical side is about a factor 2 larger than the theoretical side. Just for a rough 
estimation, let us multiplying the e\ value by a factor 2. We find that the T value at which 
the last term is 1/3 of the first term in the sum rule before differention with respect to 
1/T move to lower values 0.6 GeV and 1.9 GeV for the sum rule (|T7|) and (^) respectively. 
That is, the convergence becomes faster. The values of Q and Q 1 change to 

G = 3.6 ±1.0 ±0.7 GeV" 2 , 1.0 GeV < T < 1.5 GeV , 

G' = 1.2 ±0.1 ±0.5 , 1.8 GeV<T<2.1 GeV. (35) 

On the other hand the sum rule (El) is not good in stability. Especially, Q 1 changes 
sign near T ~ 1.2 GeV. Therefore, we do not use it in the numerical calculation of Q' . 
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The decay widths in the leading order are 

r(£* - Bn) = ^-Q 2 \q\ 5 , 

2U7T 

r(s* - sv) = ^ 2 M 5 , 
r(^ - btt) = ^-g' 2 \q\ , 

T(B[^B*n) = ^-g' 2 \q\ J (36) 

where sum over charged and neutral pion final states has been included in the above 
formulas. 

In the following we apply the leading order formulas obtained above to the excited states 
of charmed mesons. The value q for the first three processes in fl3"6|) are calculated from the 
experimental mass values of the relevent particles [JTTJ] and that for the last two processes 



are calculated from the difference between the theoretical mass values of the two doublets 
in the leading order obtained in |J and ||. For the excited states of D we obtain the 
results listed in table |, where the numbers in the parentheses are results for e\ = —0.028. 



The experimental data [[LI] and the results of || obtained with relativistic B-S equation are 
also included. For the D series the 0{1/w,q) corrections are not negligible. Yet, assuming 
the total width of D^ is dominated by one pion decay, the calculation result for D\ is in 
reasonable agreement with experimental data. The calculation result for Y[D\ — ► D*n) in 
the leading order is significantly smaller than the experimental data for T(Di — *> all). This 
can partly be attributed to mixing of the two 1 + states in the order of 0(1 /tuq) 0. This 
effect can be calculated within our formalism by introducing the 0(1/7tiq) terms in the 
Largrangian of HQET. In reference [Q], it was estimated that the T(Di — > Dp + D*p) ~ 6 
MeV, r(£>2 — > Dp + D*p) ~ 3 MeV. Therefore, the effect of mixing of two 1 + states in the 
order O^I/tuq) needs not to be very large. The calculated results for the decay widths of 
D' Q and D[ in the leading order is too small. This may due to the slow convergence of Eq. 
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Appendix: The OPE of matrix element (7r i (q)\T(q(x)q(0))\0) 

(7^)1^(0)^(0)10) = <(/ l75 + U lh U (|) afe , (37) 
(7r l (g)| J D M g«(0)g5(0)|0) = {( 9l + g 2 4)q^ 5 + h lllllh + h 2 ia^ 5 q»} aP , (38) 

(^(3)1^^(0)^(0)10) = i{[g^{ ai + a 2 i) - ia^h + ie^Yq^h 

+ {q»lv + Qul^ci + (q^ u ~ q v l»)di 
+ i(q^\v + <?^a m )<7 A c 2 + i{q^\v - q l yCr Xfl )q x d 2 
+ (ei + e 2 i)q^q v ]^} aP {^j ^ ; (39) 

fi = ~(QQ) , /a = -y; ( 4 °) 

1 i 

9i = ^(qq)o, 92 = ^, h = 0, h 2 = -—^-( qq ) o; (41) 



ai = --^T(qq)o , a 2 = —^fnml, c 1 = ^-f 7r m 2 1 , c 2 = --^—(qq) , (42) 
16/tt 36 36 24f w 

where (vr l (g)|g £> 2 7p7 5 g|0) = -if n m\ q p . 

e 1 (n-q) 2 = l -{n i (q)\q^(n-D) 2 Kq\0) , for n 2 = , (43) 
e 2 (n.q) 3 = l -(n*(q)\q^ ^i(n-D) 2 q\0), for n 2 = , (44) 
h = y , &2 = -3ci , di = , d 2 = -y + ~^(qq)o ■ (45) 

Here / w = 93 MeV. 
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Figure Captions 

Fig. 1. Dependence of Q and Q' on the Borel parameter T. The three 
curves correspond, from top to bottom, to the sum rules flllij), (|23|) and 
(|24|), respectively. 
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TABLES 

TABLE I. The decay widths of D±, D 2 , D[ and D' . The numbers in the parenthesises are 
results obtained by using e\ = —0.028 instead of —0.014. 



Transition 


Calculated Width(MeV) 


Observed Width(MeV) 


Result from || 


T(D l - 


-> D*tt) 


7.0(2.8) 




8.0 


T(D l - 


-> all) 


18.91SJ 


T(D* 2 - 


-> Dtt) 


18.2(7.3) 




15 


r(D|- 


-> D*tt) 


6.9(2.8) 




6.8 


T(D* 2 - 


-> all) 


23±| 


r(D' - 

T(D[ - 


-> Dtt) 
-> L>*tt) 


45(64) 
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